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I. INTRODUCTION 



Nine real dimensionless parameters characterize nonbirefringent modified Maxwell theory 
minimally coupled to the standard D.rae theory of massive spin-1 particles M. Recently, 
new direct laboratory bounds on these parameters have been published [4H6|], specifically, 
upper bounds on the absolute values of the isolated nine parameters. The direct laboratory 
bounds for the five parity-even anisotropic parameters are now only 1 order of magnitude 
above the previous indirect bounds [7, 8]. The direct laboratory bounds for the three parity- 
odd anisotropic parameters are still 5 orders of magnitude above the indirect bounds and the 
situation is even more dramatic for the single isotropic parameter Q. (The qualifications 
'direct' and 'indirect' will be explained later.) 

With the laboratory bounds steadily improving , we are, therefore, motivated to 

see what the chances are to tighten our indirect bounds in the coming decades. In partic- 
ular, we will focus on the present and future capabilities of the Pierre Auger Observatory 



(PAO) 111], the High Energy Stereoscopic System (HESS) of imagin g at mospheric Cherenkov 



telescopes M , and the Cherenkov Telescope Array (CTA) project U As regards g amma- 
ray telescopes, our focus is on HESS and CTA, rather than MAGIC 14| and VERITAS 15], 
because the first two offer the best prospects for detecting photons with the highest possible 
energies, which will turn out to be the crucial ingredient for one of our potential bounds. In 
that respect, also future extensive-air-shower arrays look promising fit!]. 1 

Returning to modified Maxwell theory, it needs to be mentioned that, in principle, 
the dimensionless Lorentz-violating parameters of the photon sector can be of order 1 if 
they arise from a nontrivial small-scale structure of spacetime |23[|. The single isotropic 
Lorentz-violating parameter /t tr , in particular, has a natural interpretation, if positive, as 
the excluded-spacetime- volume fraction of "defects" embedded in flat Minkowski spacetime. 
It is, thus, of fundamental importance to get as tight bounds as possible, preferably for all 
nine dimensionless Lorentz-violating parameters of the nonbirefringent photon sector but 
certainly for the isotropic one. 
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3, 



which are 



1 An entirely different class of "cosmic messengers" consists of high-energy neutrinos 
particularly well suited for the search of new Lorentz-violating and possibly CPT-breaking effects that lie 
outside the mass sector 20-221. 
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II. THRESHOLD CONDITIONS 



An interesting suggestion for obtaining bounds on the Lorentz- violating (LV) parameters 



of a hypothetical nonrelativistic theory is the following 24H26| : 



(a) with modified dispersion relations, new decay channels may appear which are absent 
in the standard relativistic theory; 

(b) this leads to rapid energy loss of particles with energies above threshold [the threshold 
energy E t ^ depends on the LV parameters k of the theory, typically with E^k) —> oo 
for k — y 0]; 

(c) observing these particles implies that they necessarily have energies at or below thresh- 
old (E < £th)j which, in turn, gives bounds on the LV parameters of the theory. 

This suggestion relies, of course, on the proper calculation of the threshold energies and decay 
rates. See, e.g., Re 
relations and Ref. 



? . |27l| for a general discussion of physically admissible modified dispersion 
28j for a general discussion of Lorentz- violating decay processes, starting 
from simple scalar models. 

Consider, now, the photon (7) of nonbirefringent modified Maxwell theory [J 2)] minimally 
coupled to standard massive spin-| Dirac particles (e.g., the electron e and proton p) in a 
particular preferred frame, for definiteness, identified with the sun-centered celestial equa- 
torial frame. 2 Two Lorentz-violating decay processes, in particular, have been considered: 



vacuum- Cherenkov radiation (p + — > p + ~f) [9 



29 



and photon decay (7 — > e + e ) 



23. 



We start by discussing the first process, vacuum- Cherenkov radiation of a proton p or heavy 
nucleus N, both, in first approximation, considered as charged pointlike Dirac particles. 

Assume the detection of an ultra-high-energy-cosmic-ray (UHECR) primary with energy 
E pTim , mass M prim , and flight-direction unit vector q prim . In terms of the standard-model- 
extension parameters [sj], the vacuum- Cherenkov process (N 
following threshold condition 0, 0, T 



pnm -> N prim 7) then gives the 



30j on the LV parameters of nonbirefringent modified 



2 More precisely, the theory considered is the one of the standard model with the experimentally determined 
values for the coupling constants and Higgs vacuum expectation value, to which a single nonbirefringent 
Lorentz-violating kinetic term of the U(l) gauge bosons is added, namely, the ks term of Eq. (16) in 
Ref. or the k' 1 ) term of Eq. (C2a) in Ref. [9]. As remarked in Sec. V of Ref. 0, the bounds discussed 
in that article and this one may also apply to theories with additional Lorentz violation in the quark and 
lepton sectors, but, then, a precise formulation of the bounds becomes rather involved. For this reason, 
the discussion here is restricted to the theory with a single Lorentz-violating kinetic term of the Abelian 
gauge field. 
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Maxwell theory: 

R[2 n tr — gp rim (k q+ ) 6c e abc — gp rim {^ c -) ab Q'prim] < {M pr i m c 2 /E prim ) , (1) 

where each repeated Cartesian index a, b, and c is summed over 1 to 3, the quantity e abc 
stands for the totally antisymmetric Levi-Civita symbol with normalization e 123 = 1, and 
the ramp function R is defined by R[x] = (x + \x\ )/2. The velocity c = 299 792 458 m/s 
on the right-hand side of ([I]) corresponds, in the Lorentz-violating theory considered, to the 
maximal attainable velocity of the Dirac particles. The nine LV parameters are collected in 
the antisymmetric 3x3 matrix k 0+ (three parameters) , the symmetric traceless 3x3 matrix 
k c _ (five parameters), and the number K tr (one parameter). 

The minimal set of UHECR events needed to obtain significant bounds from ([1]) on all 
nine LV parameters has six events (labeled I — 1, . . . , 6) with approximately equal leverage 
factors (Ei/Mic 2 ) 2 3> 1 and approximately (anti)orthogonal flight directions q ; . In an 
appropriate Cartesian coordinate system, the ideal flight directions are given by q 12 = 
(±1, 0, 0), q 34 = (0, ±1, 0), and q 56 = (0, 0, ±1). Strictly speaking, three orthogonal 
events suffice for the parity-even parameters K e _ and a single event for the positive isotropic 
parameter K tr . 3 As emphasized in Ref. 25] and reiterated in Refs. jsiH], these bounds require 



as only input the mere existence, at the top of the Earth's atmosphere, of charged cosmic-ray 
primaries with travel lengths of a meter or more. Hence, these bounds are independent of 
the distance to the (astronomical) source. 

Next, restrict to the isotropic sector of nonbirefringent modified Maxwell theory and 
consider the case of a small negative parameter K tr . Then, vacuum- Cherenkov radiation 
is no longer possible, but another decay process becomes available, namely, photon decay 
into an electron-positron pair. The photon-decay process (7 — > e + e~) gives the following 
threshold condition 0, Q, 25] on the isotropic LV parameter of modified Maxwell theory: 



R[ -«tr] < 2 (M e c 2 /£ 7 ) 2 , (2) 

with E 1 the energy of the primary photon, M e the mass of the electron, and R the ramp 
function defined below ([[]). 

The minimal set of gamma-ray events needed to obtain a significant bound from ([2]) 
contains a single event with leverage factor {E^jM e c 2 ) 2 ^> 1 and an arbitrary flight direction 



3 Remark that, with so few events, the bounds are somewhat weakened by the appearance of geometrical 
factors in the function argument on the left-hand side of (TT]). The geometrical factor of the second (parity- 
odd) term, for example, is of order cos7r/4 = l/-\/2 ss 0.71. See also the penultimate paragraph in Sec. II 
of Ref. 7] for further discussion. 
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(the considered Lorentz violation being isotropic). Again, the only input for this bound is 
the mere existence, at the top of the Earth's atmosphere, of a primary photon with a 
travel length of a meter or more, making the bound independent of the distance to the 
(astronomical) source. 4 



III. PROTOTYPICAL DATA SAMPLES 



Two fiducial data samples are presented in this section, the first with UHECR events 
for application of the vacuum- Cherenkov threshold condition ([1]) and the second with TeV 
gamma rays for application of the photon-decay threshold condition (J2]). Our aim is to be as 
concrete as possible, but it is clear that the data samples considered are only examples and 
can certainly be modified or improved in the future, as will be discussed further in Sec. IV Al 



For the first prototypical data sample, we consider hybrid UHECR events 3JJ detected 
by Auger (short for the Pierre Auger Observatory, PAO). Recall that a so-called 'hybrid' 
Auger event refers to an extended air shower which has been observed simultaneously by 
the two typ es of detectors, the water- Cherenkov surface detectors and the fluorescence tele- 



scopes 



Specifically, the following sample £jv' of N' = N/2 hybrid Auger events is obtained from 



the two energy bins below E = 10 19 eV in Figs. 2 and 3 of Ref. [3l| : 

N/2 = 131 + 96 = 227, (3a) 

E ~ 8 EeV , (3b) 

< A max > ~ 750 g/cm 2 , (3c) 

rms(A max ) ~ 50 g/cm 2 , (3d) 

where the reason for writing N/2 will become clear shortly and where < X max > stands for 
the average of X max , the atmospheric depth of the shower maximum (here, determined by 
the fluorescence telescopes of Auger). The relative uncertainty of the energy measurement 



3b]) is taken equal to 25% (slightly larger than the value of 22% quoted in Ref. [31]). 



4 If it could be established that certain specific astronomical sources emit UHECR protons or multi-TeV 
photons which definitely radiate their energy away or decay before they reach the Earth, then a finite 
range of nonzero values of the relevant LV parameters could be obtained, assuming the LV effects to be 
responsible for the observed disappearance of the particles. Such a range could be expected to have a 
(mild) dependence on the source distances, whereas the bounds of the present article are strictly source- 
distance independent. 
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Remark that the fluorescence telescopes of Auger [ll|, [3lj] partly track the primary particle 
and its secondaries along their way in the atmosphere, following them over a length of several 
kilometers. 

The sample Ejv/2 from (J5a} -( l3d]) is, however, not quite isotropic, as the celestial coverage 
of Auger-South is not complete. Awaiting results from Auger-North (or a similar cosmic-ray 
observatory on the northern hemisphere), we can proceed as follows. An artificial sample 
S^r is constructed by adding to each event of T, N / 2 one in the opposite direction: 

S w = { {(E n , q n , C n ), (E n , -q n , C n )} I n = 1, . . . , N/2 } , (3e) 

where q n is the flight-direction vector of event n and C n contains further characteristics such 
as A" max n . In this way, the sample T, N covers the whole celestial sphere, even though certain 
regions may be covered somewhat more densely than others. [There are, of course, other 
ways to make the sample £at/2 more isotropic, for example, by keeping the same number of 
events but randomly flipping the flight-direction vectors.] 

The last two assumptions are that the sample Ejy contains N p 3> 1 events with proton 
primaries and that the arrival directions of these N p events are representative of an isotropic 
distribution. Specifically, we assume that 

N p > 48, (3f) 

and that the proton arrival directions cover the celestial sphere completely and more or 
less homogeneously. The isotropy condition can be described in words as follows: to have 
approximately N p /2 events in each hemisphere (solid angle 2tt sr), approximately N p /4 
events in each quadrant (solid angle 7r sr), approximately N p /8 events in each octant (solid 
angle 7r/2 sr), etc. Restricting to octants, one possible mathematical formulation of the 
isotropy condition for the proton arrival directions is 

V n e{i,2, ...,8} : jVp(An n = 7r/2) > 1, (3g) 

which must hold for any choice of octants AQ n that cover the celestial sphere (Q = Att) 
completely. Condition fl3g"D is not too difficult to satisfy with the relatively large total number 
of protons from (|3fj). For the data sample T, N , conditions (J3f|) and fl3g| ) could perhaps result 
from those events with A max > 825 g/cm 2 , according to the calculations shown in Fig. 3 of 
Ref. Q (see also Table I of Ref. Q). 

The basic idea, now, for obtaining new bounds from ([T]) and the prototypical data sample 
(J3]) is that, without knowing precisely which of the iV = 454 events would correspond to 
a proton, the knowledge would suffice that there would be a subset with N p 1 proton 
events distributed isotropically. Clearly, this requires the reliable determination of the proton 
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fraction P p in the sample considered (with the certainty that P p > 0) and the information 
that the protons in this UHECR sample have an overall isotropy of the arrival directions 
(irrespective of possible small-scale clustering). In short, the challenge for Auger will be to 
identify a 10-EeV UHECR sample with enough isotropic protons in it. Further discussion 
will be given in Sec. IV Al 

Turning to the second prototypical data sample, recall that 80 TeV gamma rays have 



been detected by HEGRA from the Crab Nebula 32] and 10 2 TeV gamma rays by HESS 



from the shell-type supernova remnant RX J1713. 7-3946 [33j, |34|. The relative uncerta inty 



of the energy measurement for this last extended source equals 20 % according to Ref. [33]. 



Remark that the Cherenkov light detected on the ground [35] indirectly tracks the primary 
photon and its secondaries along their way in the atmosphere, essentially following them 
over a length of several kilometers. 

Specifically, the second sample is taken to have a single fiducial photon event with an 
arbitrary arrival direction, 

N 7 = 1, (4a) 
£ 7 ~ 10 2 TeV. (4b) 

Here, the crucial assumption is that at least one such photon has been detected unambigu- 
ously (or, more realistically, that this photon is one of a larger sample of photon candidates 
identified with high probability). 



IV. POTENTIAL SENSITIVITIES 

Relying on the proton primaries contained in the prototypical UHECR sample ([3]), the 
vacuum- Cherenkov threshold condition ([1]) could be used for the numerical values 

E pvim = (8 ± 2) EeV , M prim = M p = 938 MeV , (5) 

explicitly showing the one-sigma error in the energy determination and inserting question 
marks to emphasize that the data sample is still artificial. Assuming a sufficiently large 
number N p of proton primaries, the directions q pnm) ; could be taken optimal, so that the 
geometrical factors entering the argument on the left-hand side of ([I]) would be close to 
their maximal values. The main contribution to the statistical error on the bounds of the LV 
parameters would come from the energy uncertainty and would correspond to approximately 
50 % for this sample. 
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The resulting two-cr indirect bounds on the nine isolated LV parameters of nonbirefringent 
modified Maxwell theory would be: 

\(Ko + ) (ab) \ (ahMmY(12 , < 1-4 xl(T 20 , (6a) 



(06) = (23), (31), (12) 
(^c-) ( ) I (06) = (11), (12), (13), (22), (23) < 3 X 10 °' 



K tT < 1.4 x 10" 20 , (6c) 

which, for the first two entries, would improve by 2 orders of magnitude and, for the last 
entry, by 1 order of magnitude upon the previous bounds in Ref. js]. There would be no 
mystery where the improvement would come from: comparing the values (jSJ) to those used 
in Ref. js], the gain from a much smaller mass value M prim would overcome the loss from a 
somewhat larger value of the inverse energy 1/E pr i m . 

With the prototypical gamma-ray sample fll]), the photon-decay threshold condition (J2J) 
could be used with numerical values 

E 1 = (100 ± 20) TeV, M e = 0.511 MeV. (7) 

The resulting two-cr indirect lower bound on the single isotropic LV parameter of non- 
birefringent modified Maxwell theory would be 

n tr > -0.9 x 10~ 16 , (8) 

which would improve by 1 order of magnitude upon the previous two-cr lower bound K tT > 
—0.9 x 10 -15 [9j. 5 Again, there would be no mystery where the improvement would come 
from, namely, the somewhat larger energy value E 1 in ([7]) compared to the one used in 
Ref. Q. 

Both indirect bounds ([6]) and (jSJ) would scale with the inverse of the energy square of the 
incoming particle, that is, with [(8 EeV)/-E pr i m ] 2 and [(10 2 TeV)/E 1 \ 2 , respectively. But, 
for the moment, these bounds are only potential bounds as they rely on certain assumptions 
(e.g., N p ^> 1) or on low-statistics data (e.g., E 1 = 10 2 TeV). 



V. DISCUSSION 

A. Experimental issues 

As mentioned in the previous section, bounds fl6]) and jSJ) are, for the moment, only 
indicative of the potential sensitivity of Auger and HESS/CTA, because they either rely on 

5 The qualitative lower bound on K t r at the — 10~ 16 level (no confidence level specified) from Ref. J36| relies 
on a detection of 50 TeV gamma rays from the Crab Nebula, which is of marginal significance [37| . 
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certain assumptions or use low-statistics data. We expect that both bounds will be realized 
in the coming decades, derived without additional assumptions and from high-statistics data. 
Let us now discuss the arguments on which this expectation is based. 

The potential UHECR bounds (E]) rely on two main assumptions for the event sample 
considered: (i) the large-scale isotropy of the UHECR arrival directions (applicable to all 
primary subspecies), and (ii) the presence of a nonvanishing fraction of proton primaries 
or, at least, of light-nucleus primaries. Both assumptions will be tested (and, most likely, 
verified) with further Auger-South data and forthcoming Auger-North data. 

At this point, it is already clear that it will be of critical importance for Auger to obtain, 
for a subset of events, a reliable determination of the type of primary (e.g., a hydrogen, 
helium, oxygen, or iron nucleus), preferably by use of different types of diagnostics (for 
example, in addition to A max , number densities of muons and arrival times of particles 
reaching the surface {38j]). As mentioned in Sec. [TT1 only six ( ant i) orthogonal events with 
(M p r i m c 2 /E prim) 2 ~ 10~ 20 would suffice, as long as M prim would be determined accurately. 

Alternatively, a nonzero fraction P; of primary type % (for example, % — p, He, O, or Fe) 
would suffice for establishing the same bound, provided a large enough sample of isotropic 
events would be available with (MjC 2 /Pj) 2 ~ 10~ 20 [this is, in fact, the strategy followed 
with the prototypical data sample ([3]) used to obtain the potential bounds for i = p and 
Mi = 0.938 GeV]. This strategy, with a nonzero Pi for an appropriate sample of UHECR 
events, appears the most promising for improving the bounds of Ref . [8] . 

Recent observations from Auger 3l| indicate an increasing fraction of heavy nuclei with 
increasing primary energy E. Still, it is very well possible that the fraction of light nuclei 
(e.g., hydrogen) is small but nonzero at any fixed value of E. In fact, a recent suggestion [39] 
to explain the rising fraction of ironlike nuclei seems to imply just that (e.g., P p ~ 10 % 
and Pf e ~ 90 % at E ~ 8 EeV for the simple model of Fig. 1 of Ref. [39j). It, then, remains 
to determine the optimal primary type i with the smallest value of MiC 2 /Ei and enough 
events (Pi N > 10) having arrival directions more or less isotropically distributed over the 
celestial sphere (irrespective of the positions of the original UHECR sources). As mentioned 
already in Sec. IHH the argument for obtaining the bound is statistical, without the need to 
know precisely which event of the iV-event sample corresponds to species type i, as long as 
we can be sure about the nonzero value of the fraction Pi and the arrival-direction isotropy. 

Now turn to the potential gamma-ray bound (JHJ), which relies on the detection of 10 2 TeV 
gamma rays by HESS. For the moment, this detection has only a marginal significance 
(less than two-a according to Fig. 3 and Table 5 of Ref. 33]). Future improvement can 
be expected from further HESS (or VERITAS) observations and certainly with CTA fl3|. 
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which may have a 10 times greater sensitivity than HESS at energies E 1 ~ 10 2 TeV. If 
CTA (or an extensive-air-shower array 16|) unambiguously detects 10 2 TeV photons from 
whichever astronomical source, the potential bound flH]) will be transformed into one of the 
cleanest bounds on isotropic Lorentz violation in the photon sector. 

In addition, there is always the possibility of a detection by Auger of ultra-high-energy 



photons 40j, |4l|. Note that a photon event with E 1 = 10 18 eV, for example, would improve 
the lower bound ©, based on fl2} and ©, by a factor 10 8 to the level of -0.9 x 10~ 24 . 



B. Theoretical issues 

Following up on these experimental considerations, we have two remarks on theoretical 
issues, one regarding the distinction of 'direct' vs. 'indirect' bounds already made in Sec. |T] 
and the other regarding the distinction of 'earth-based' vs. 'astrophysical' bounds. 

First, a direct bound on the LV parameters of nonbirefringent modified Maxwell theory []]- 
is taken to refer to experiments which directly test the modified propagation properties 
of the nonstandard photons, for example, by measuring positions and arrival times of pulses 
of light. An indirect bound, instead, is taken to refer to experiments which look for derived 
effects resulting from the nonstandard dispersion relations of the photons. These derived 
effects [here, vacuum- Cherenkov radiation and photon decay] are, however, unambiguously 
calculable from the theory considered in this article, which has a modified pure-photon 
sector and a standard matter sector with minimal photon-matter couplings (governed by 
the principle of Kauge invariance). In fact, the energy thresholds used for the indirect 
bounds of Refs. [8|, |9J and the present article only rely on energy-momentum conservation. 
For this reason, we do not consider indirect bounds to be necessarily less reliable than direct 
bounds. 

Second, it is important to realize that the existing indirect bounds from Refs. 

as 

and the 

potential indirect bounds from the present article are earth-based bounds (the atmosphere 
of the Earth corresponding to the 'front end' of the detector [ill Q) an d n °t astrophysical 
bounds, which may rely on additional assumptions about the astronomical source and the 
propagation distance of the primary particles. We have already mentioned in Sees. [TT] and 
IIHI that our indirect LV bounds require the existence of primaries with track lengths of a 
meter and that Auger measures shower tracks which start at a point high up in the Earth's 
atmosphere and run for several kilometers. Hence, we can simply repeat what was written in 
Ref. [8 1, "these Cherenkov bounds only depend on the measured energies and flight directions 
of the charged cosmic-ray primaries at the top of the Earth atmosphere." The same holds 
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for the photon-decay bound, now referring to measurements of the photon primary. As such, 
these earth-based bounds are not inherently inferior to laboratory bounds. 



C. Comparison of bounds 

Tabled] now compares the different types of bounds for the parameters of nonbirefringent 
modified Maxwell theory minimally coupled to standard Dirac theory. The second data 
row of this table refers to a potential sensitivity of order 10~ 20 [5] which may be achieved 



with cryogenic resonators 10]. As it stands, direct laboratory experiments are affected by 
the earth-velocity boost factor ~ 10~ 4 , which reduces the sensitivity for the parity-odd 
parameters k 9 + linearly and for the isotropic parameter K tr quadratically. There has been 



a suggestion [42| of how, in principle, this boost-dependence problem may be overcome for 
the parity-odd parameters, but, to the best of our knowledge, this has not yet been realized 
experimentally. 

As the effects of the Earth's velocity are negligible for the indirect bounds from vacuum- 
Cherenkov radiation and photon decay j^oj], the corresponding indirect bounds for the 
parity-odd and isotropic parameters are much stronger than the direct bounds. 6 Note also 
that it is entirely possible that future data samples from Auger-South+North and CTA may 
provide us with 3 times larger energies than used in (J5]) and (J7|), so that the sensitivities of the 
last data row in Table [I] would be improved by an order of magnitude to the level of 10~ 21 for 
the eight anisotropic parameters 7 and (— 10~ 17 , +10 -21 ) for the single isotropic parameter. 
Especially, the latter projected CTA/Auger sensitivities for the isotropic parameter K tT would 
be a welcome complement to future laboratory bounds. 



VI. CONCLUSION 



The current direct laboratory bounds 
nonbirefringent modified Maxwell theory 



4-6] on the nine Lorentz- violating parameters of 
l|-|3) minimally coupled to standard Dirac t 



are absolutely remarkable. Still, they fall short of what has been achieved indirectly 




The same holds for an indirect laboratory bound (relying on derived effects for the synchrotron radiation 



rate at LEP) which gives |«t r | < 5 x 10 15 at the two-cr level [43 1. 



7 Inferred two-cr bounds on the eight anisotropic parameters already exist at the 10 20 level (even better 



for the parity-even parameters) from an atomic-fountain-clock laboratory experiment 44( , if a coordinate 
transformation is used to move the isolated Lorentz violation in the proton sector (eight parameters) to 
the electromagnetic sector (eight parameters). For further discussion of these coordinate transformations, 
see, e.g., App. A in Ref. [45[ and Sec. VII in Ref. 



11 



TABLE I: Orders of magnitude for existing and potential two-c bounds on the nine Lorentz- 
violating parameters of nonbirefringent modified Maxwell theory [l-jj (minimally coupled to stan- 
dard Dirac theory) from direct laboratory and indirect earth-based experiments. The considered 
indirect bounds for the isotropic parameter /ttr have different orders of magnitude depending on 
the sign of /5 tr [the values of the last data row, for example, correspond to ([5]) and (|5"cj) in the main 
text]. The second data row for potential direct laboratory bounds shows explicitly the effects from 
the earth- velocity boost factor /3m ~ 10~ 4 . 



Type of bound 



K o+ 



Reference 



Existing, direct 
Potential, direct 
Existing, indirect 
Potential, indirect 



1(T 17 

10 -20 ? 

io- 18 

1Q -20 ? 



io- 13 

10 -20+4 ? 

io- 18 

1Q -20 ? 



io- 8 

1Q -20+8 ? 

-io- 15 , +10- 19 ) 

- 10~ 16 ?, +10" 20 ?' 



[4-6] 
[5, 1Q] 
[8,9] 
This paper 



as can be seen by comparing the first and third data rows in Table [H Instead of arguing 
about the relative merits of these direct and indirect bounds (cf. Sec. IV BP . we adopt the 
following pragmatic approach. 

The existing indirect earth-based bounds of Refs. js, 9] can be used to predict that future 
direct or indirect laboratory experiments will not detect nonbirefringent Lorentz violation 
in the photon sector at levels of order IO -18 for the eight anisotropic parameters and 10~ 15 
for the single isotropic parameter. Furthermore, this article has shown that experiments 
at present and future observatories such as Auger and CTA have the potential to reduce 
these predictions to the 10 -20 level for the anisotropic parameters and to the 10~ 16 level 
for the isotropic parameter (see the last data row in Tabled]). Perhaps the most interest- 
ing predictions of Auger and CTA would be for the isotropic parameter K tr , which might 
ultimately drop to levels of order (— 10~ 17 , +10~ 21 ) for energy values E pr i m = 25 EeV and 
E~ = 3 x 10 2 TeV replacing those of Eqs. ([S]) and ([7]), respectively. 
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